clear all;close all

A = xlsread('pitch_angle_scattering.xlsx');

figure
bar(A(:,1),A(:,2));
set(gca,'fontsize',20);
xlabel('Pitch angle [^\circ]');ylabel('Number of test electrons');
axis square
title(['\omega=2.1\omega_{ce}, s=+1, E_p(0)=39 keV, t=0.4 ms']);
axis([57 104 0 60])

figure
bar(A(:,1),A(:,3));
set(gca,'fontsize',20);
xlabel('Pitch angle [^\circ]');ylabel('Number of test electrons');
axis square
title(['\omega=2.1\omega_{ce}, s=-3, E_p(0)=1.15 MeV, t=0.4 ms']);
axis([57 104 0 240])

figure
bar(A(:,1),A(:,4));
set(gca,'fontsize',20);
xlabel('Pitch angle [^\circ]');ylabel('Number of test electrons');
axis square
title(['\omega=0.597\omega_{ce}, s=+1, E_p(0)=39 keV, t=0.4 ms']);
axis([57 104 0 140])

figure
bar(A(:,1),A(:,5));
set(gca,'fontsize',20);
xlabel('Pitch angle [^\circ]');ylabel('Number of test electrons');
axis square
title(['\omega=0.126\omega_{ce}, s=+1, E_p(0)=1.15 MeV, t=0.4 ms']);
axis([57 104 0 420])